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FOREWORD 
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The inflight short wavelength (<2000 m) gust environment is critical to 
the prelaunch simulation and flight evaluation of a space vehicle. It is 
desirable to have an accurate estimate of environmental wind at launch time. 
However, due to the limitation of existing operational wind prediction tech- 
niques and constraints imposed by detailed wind profile measurement and data 
processing procedures, it is impossible to state the inflight wind environment 
that will occur at launch or prior to launch with absolute certainty. Accord- 
ingly, the risk of exceeding critical inflight wind conditions and thus 
critical vehicle responses must be assessed with statistical procedures. 

In this report the statistics of inflight gustspectra will be examined 
in order to define the changes that can occur in the short wavelength gust 
environment. Such statistics must be determined empirically and this is the 
philosophy used in the analysis described herein. The short wavelength gust 
spectra used in this study were derived from FPS-16 Radar /Jimsphere detailed 
wind profiles. The data consistad of approximately 2700 wind profiles measured 
at Cape Kennedy, Florid? during the period 28 November 1964 to 11 May 1967. 

The smallest time period between the various profiles was approximately 
3 hours. Thus, it was possible to develop wind spectra time change statistics 
for time periods as short as 3 hours and as long as desired. The largest time 
period used in this study was 72 hours. The wind gust profile spectra were 
calculated with standard correlation-Fourier transform techniques as described 
in reterence 1. The inflight wind gust profiles were obtained by removing the 
long wavelength (> 2000 m) Fourier components. This was accomplished using 
the high-pass digital filters described by DeMandel and Krivo (ref. 2) which 
were primarily developed to filter geophysical data of the kind under con- 
sideration* 

The time change statistics of the inflight gust spectra developed give 
valuable information about the persistency of gust amplitude in time. In 
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general, the persistence of the various Fourier components of the wind field 
tends to be a monotonically increasing function of horizontal wavelength, i.e., 
the longer the wavelength of the wind feature the more persistent it is in 
time. For example, synoptic scale cyclones with horizontal wavelength on the 
order of 2000 km tend to have a typical persistence in a region of 3 to 5 days, 
while a small turbulent eddy, with a horizontal wavelength in the order of 100 
meters, for example, has a characteristic time scale or persistency on the 
order of a minute. However, the variations or wavelengths considered in this 
report are along the vertical rather than the horizontal as in the particular 
examples noted above. The largest variation or wavelength along the vertical 
is on the order of the height of the tropopause ('10-15 km) which should be 
contrasted with the typical horizontal synoptic scale of 2000 km noted above. 
Nevertheless, the monotonically increasing relationship between wavelength and 
persistence as discussed in the context of horizontal variations or wavelengths 
is still valid in a general sense for the vertical wavelength Fourier compon- 
ents of the wind field. This results from the fact that as the horizontal 
scale of atmospheric flows increase so does the associated vertical scale. 

Thus, for example, the typical vertical extent of synoptic scale eddies (hori- 
zontal wavelength '1000-3000 km) is typically on the order of 10-15 km, while 
the typical vertical extent or waveler.gth of the mesoscale (horizontal wave- 
length '10-1000 km) motions of the atmosphere is on the order of 1-5 km and 
the associated microscale (horizontal wavelength < 10 km) vertical wavelengths 
are on the order of 100-1000 meters. 
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Section II 
THE DATA 

The analysis herein is based on a set of statistics of time- •: of 

detailed wind profile spectra provided by the Aerospace Environment Divi.uen of 
the NASA-Marshall Space Flight Center. These spectra were calculated with 
detailed wind profile data acquired with the l’?S-16 Radar/Jimsphere system at 
Cape Kennedy, Florida during the period of 28 November 1964 to 11 May 1967. 

A discussion of this wind sensing technique can be found in references 3 and 
4. The empirical spectral time change statistics provided by NASA Included 
distribution functions at various percentage points for the quantity 


C(K;t,-r) - $(K;t + x) - <fr(lC;t) 


( 2 - 1 ) 




m 


where (K; t ) is the spectrum at vertical wave number K of the vertical wind 
profile of horizontal wind speed at time t, and t is the time Interval between 
wind profiles. The elapsed time or time lag t took a value between 3 to 72 
hours (3, 6, 7.5, 9, 12, 18, 24, 48, 72 hours) and the quantity K range between 
0.00025 cycles m 1 to 0.0075 cycles m * (0.00025, 0.0005, 0.001, 0.0015, 0.002, 
0.0025, 0.00375, 0.005, 0.00625, 0.0075 cycles m ^). In addition, the empiri- 
cal statistics included associated unbiased values of the mean, variance, 
skewness and kurtosis calculated according to the following formulae 


mean 


1 "i ’ i £ 


1/2 

standard deviation: o * p ' 




skewness: S » — 


*3 (n-1) (n-2) U i ’ y i r 


( 2 - 2 ) 


(2-3) 


(2-4) 


2-1 




i 

h. 
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kurtosis: K 


1_ f n 2 - 2n - 3 r 
o 4 \(n-l) (n-2) (n-3) ^ 



4 


- 3(2n / e . 

n(n-l) (n-2) (n-3) 



For details concerning these formulae see reference 5. 

To compute the distribution function and the above moments it was assumed 
that the process £(K;t,T) is stationary and that the ergodlc hypothesis is 
available. Thus, the empirical distributions and moments of £ are functions 
of K and t only. 

In an attempt to graduate the spectral statistics the data in the sample 
were categorized according to the largest mean flow wind speed in the 5-15 km 
region at time t, which will be denoted by V. The mean flow is defined to be 
the flow associated with Foruier component possessing vertical wavelengths 
>2000 m. The categorization scheme consisted of the partition of the sample 
of values for £(K;!-,t) into two subsamples, i.e., a sample of £(K;t,i) for 
V <_ 45 msec * and one for V > 45 msec”*. Thus the distribution function and 
moments in this report are conditional ones, with the conditionalized parameter 
being the quantity V. 
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Section III 

DATA ANALYSIS 

In this section the data analysis methods that were used to construct a 
model of the time difference statistics of detailed inflight wind profile 
spectra will be discussed. Preliminary analysis revealed that it was possible 
to develop a "universal" statistical distribution of the spectral changes which 
encompassed all vertical wave numbers, wind speed categories, and more impor- 
tantly, elapsed time. This, as will be shown later, Is a fortuitous result 
which will be extremely useful in operational applications. 

To develop a "universal" probability distribution of wind spectra time 
changes the distributions of £ supplied by the NASA were transformed to dis- 
tributions of the standardized variable 


€< WV ~ VvV 


where 


p (P i ,t j ,K k ) - o(T j ,K k ) 


p is the standardized wind spectrum difference 
£ is the original wind spectrum difference 
is the mean of £ 

a is the scandard deviation of £. 


(3-1) 


; \ 


The notation Inside the parentheses serve to indicate percentage point 
(p^ « 1, 5, 10, 25, 50, 75, 90, 95, 99, 99.9 percent), elapsed time (t^ » 3, 

6, 7.5, 9, 12, 18, 24, 48, 72 hours), and wave numbers (K^ • 0.00025, 0.0005, 
0.001, 0.0015, 0.002, 0.0025, 0.00375, 0.005, 0.00625, 0.0075 cycles m* 1 ) . 

After reviewing the plot of the standardized wind spectra difference p 
versus percentage point p^ for various elapsed time tj, they showed no signif- 
icant difference. Therefore, they can be averaged into one curve by the 
following formula. 
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p(p i’V 


S e^yV « «<y V 


(3-2) 


l "(T ,K ) 
j-1 J 


where p is wind spectra difference averaged over elapsed time, n is the number 
of distributions for specified elapsed time and wave number. The notation 
inside the parentheses is previously defined. Also, plots of the wind spectra 
difference averaged over elapsed time for various wave numbers show that the 
probability distribution curves of the wind spectra difference averaged over 
elapsed time did not vary much among various wave numbers. A similar formula 
can be used to average these curves to obtain the "universal" probability dis- 
tribution of wind spectra time changes. 

10 

I x 

(3-3) 


Pfpp 


k=l 


10 

l * 0 ^) 


k=l 


Al. 

\ 


where 


p(p^) is wind spectra difference averaged over wave number and elapsed 
time 

ra(K^) is the number of distributions for wave number and 
9 

m(Kk> * j^l n<T j' K k ) ‘ 


\ 

i 


y, 

* 




Thus, the "universal" probability distribution of wind spectrum change 
has been averaged in terms of the random variable p. To recover the K distri- 
bution, models are required for and o. The mean value, u^, is identically 
equal zero (as will be shown later) . To develop a model of the standard devia- 
tion as a function of wave number K, the following equation was used to average 
the standard deviation of wind spectrum differences over elapsed time. The 
standard deviation o is an exponential function of wave number K, as will be 
shown later. 


3-2 
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0( V 


[ 

Jzl 


°( T j> K k ) x nCTj.K^) 


I n(T,X) 
j«l 3 * 


(3-4) 


where o(K^) is standard deviation of wind spectra average over elapsed time. 
All other notations have been previously defined. 


Also, similar formula will be used to get p^K^) the standardized mean of 
wind spectra difference averaged over elapsed time t. , 


3-3 
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Figures 4-1 through 4-10 show the probability curves of the standardized 
gust spectra difference averaged over elapsed time t for the wave number range 
0.00025 cycles m ^ through 0.0075 cycles m ^ for the wind speeds less than 
45 msec ^ . These curves match very well from the 10 percent level up to 95 
percent level, but diverge on both ends. There are no definitive patterns for 
the divergence so that it is tentatively concluded that no trends exist in the 
tails. This is valid for engineering purposes. An inspection of Figures 4-1 
through 4-10 shows that the transformation of the random variables f; to p 
eliminates the dependence on wave number. To develop a "universal" distribution 
curve of p, the curves in Figures 4-1 through 4-10 were averaged to obtain a 
composite. Table 4-1 and Figure 4-11 shows this composite averaged probability 
curve of the standardized wind spectra difference ("universal" probability 
curve) for the wind speeds less than 45 msec Figures 4-12 through 4-21 show 
the probability curves of the standardized gust spectra difference averaged over 
elapsed time T for the wave number range 0.00025 cycles m ^ through 0.0075 
cycles m ^ for the wind speeds greater than 45 msec Similar to those curves 
of the previous case these curves agree very well from the 20 percent level up 
to 95 percent level, but diverge on both ends. It is also tentatively concluded 
that no trends exist in the tails. The composite averaged probability curve of 
the standardized wind spectra difference ("universal" probability curve) for the 
wind speeds greater than 45 msec ^ is shown in Table 4-2 and Figure 4-22. 
Comparison of Figures 4-11 and 4-22 shows that the two curves agree very well 
from the one percent level up to the 99 percent level. There is a slight 
difference at the 99.9 percent level, however, this is unimportant because 
most applications will be concerned within probability levels less than or 
equal to 99 percent. 

Figures 4-23 and 4-24 show the standard deviation of the wind spectra 
difference versus elapsed times for the wave number range 0.00025 cycles m ^ 
through 0.0075 cycles m ^ for the wind speeds less than 45 msec ^ and greater 
than 45 msec * respectively. Similarly, Figures 4-25 and 4-26 show the stan- 
dardized mean of the wind spectra difference versus elapsed times for the wave 
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Table 4-1. PROBABILITY DISTRIBUTION OF STANDARDIZED WIND SPECTRA 
DIFFERENCE AVERAGED OVER ELAPSED TIMES AND WAVE NUMBER 
FOR WIND SPEED LESS THAN 45 M SEC-1 ("UNIVERSAL" 
PROBABILITY DISTRIBUTION) 


PERCENT LEVEL 

WIND SPECTRA DIFFERENCE 

1 

-1.334 

5 

-1.199 

10 

-1.050 

25 

-0.740 

50 

-0.202 

75 

0.522 

90 

1.331 

95 

1.925 

99 

3.459 

99.9 

4.919 


Table 4-2. PROBABILITY DISTRIBUTION OF STANDARDIZED WIND SPECTRA 
DIFFERENCE AVERAGED OVER ELAPSED TIMES AND WAVE NUMBER 
FOR WIND SPEEDS GREATER THAN 45 M SEC -1 ("UNIVERSAL" 
PROBABILITY DISTRIBUTION) 


PERCENT LEVEL 

WIND SPECTRA DIFFERENCE 

1 

-1.382 

5 

-1.231 

10 

-1.072 

25 

-0.701 

50 

-0.160 

75 

0.521 

90 

1.370 

95 

2.037 

99 

3.533 

99.9 

4.343 
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number range 0.00025 cycles ra * through 0.0075 cycles m ^ for the wind speeds 
less than 45 msec * and greater than 45 msec \ respectively. There are no 
definite trends in Figures 4-23 and 4-24. But, they show that the standard 
deviation of the wind spectra difference changes rapidly with respect to 
elapsed time at elapsed time less than 24 hours for all wave numbers. Also, 
in Figures 4-25 and 4-26, the standardized mean of the spectral changes do not 
have patterns along various wave numbers. They do show rapid fluctuation with 
respect to elapsed times at elapsed time less than 24 hours for all wave numbers. 
This variability is due to statistical instability. Tne mean, by necessity, 
must be identically equal to zero. To prove this, note that £ = 4> ( t + t) - <}>(t) 
ensemble average, yields 


= ^(t + T)> - <^(t)^> 


where<Q ^is the ensemble average operation. In a stationary process 
<^(t + 1)^> = so that^£^ K 0. Furthermore, ergordicity demands that 

ensemble averages equal time averages. These are precisely the assumption used 
in the analysis, so by necessity ■ 0. Tables 4-3, 4-4, and Figure 4-27 show 
the standardized mean of the wind spectra difference averaged over elapsed times. 
There are similar patterns in both cases. As shown, theoretically, the averaged 
normalized mean of the wind spectra difference should equal to zero. Thus, the 
non-zero mean values of £ in these figures represent a bias of the data. Tables 
4-3, 4-4, and Figure 4-28 show the standard deviation of the wind spectra dif- 
ference averaged over elapsed times. A power law of the form 



(4-1) 


with q = -2.567 fit the results quite well. This is shown in Figure 4-28 
where 


o(K) is the standard deviation of the wind spectra difference at wave 
number K. 

o(K ) = 0.178, is the standard deviation of the wind spectra difference 
° at wave number k “ 0.001 cycles 
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Table 4-3. STANDARD DEVIATION AND STANDARDIZED MEAN OF THE WIND 
SPECTRA DIFFERENCE AVERAGED OVER ELAPSED TIMES 
(WIND SPEEDS LESS THAN 45 M SEC' 1 ) 


K 

STANDARD DEVIATION 

MEAN/STANDARD DEVIATION 

0.00025 

5.164204 

0.035962 

0.0005 

1 .006005 

-0.396706 

0.001 

0.19115 

-1.048616 

0.0015 

0.072866 

-0.417769 

0.002 

0.019199 

-0.065735 

0.0025 

0.009953 

-0.541371 

0.00375 

0.005549 

-0.591095 

0.005 

0.003138 

0.355242 

0.00625 

0.001663 

-0.096475 

0.0075 

0.001 

-0.068066 


Table 4-4. STANDARD DEVIATION AND STANDARDIZED MEAN OF THE WIND 
SPECTRA DIFFERENCE AVERAGED OVER ELAPSED TIMES 
(WIND SPEEDS GREATER THAN 45 M SEC' 1 ) 


K 

STANDARD DEVIATION 

MEAN/STANDARD DEVIATION 

0.00025 

5.908361 

0.427843 

0.0005 

0.922322 

-0.346052 

0.001 

0.146474 

-1.503471 

0.0015 

0.055182 

-0.783911 

0.002 

0.018193 

0.008803 

0.0025 

0.009121 

-0.455056 

0.00375 

0.003837 

-1.302469 

0.005 

0.003306 

0.645182 

0.00625 

0.002606 

0.503403 

0.0075 

0.001302 

0.339731 
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It should be noted that Figures 4-11 and 4-22 show negative values for j 

57 percent of the occasions. This means the wind spectra decreases 57 percent ! 

I 

of the time, and thus it might be concluded that the Eulerian increases in j 

gust spectra are faster than the decreases. < 

\ 
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Figure 4-2. PROBABILITY DISTRIBUTION CURVE OF STANDARDIZED WIND SPECTRA 
DIFFERENCE AVERAGED OVER ELAPSED TIMES AT WAVE NUMBER 0.0005 
CYCLES M-l FOR WIND SPEEDS LESS THAN 45 M SEC-1 
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PROBABILITY (%) 


Figure 4-3. PROBABILITY DISTRIBUTION CURVE OF STANDARDIZED WIND SPECTRA 
DIFFERENCE AVERAGED OVER ELAPSED TIMES AT WAVE NUMBER 0.001 
CYCLES M“i FOR WIND SPEEDS LESS THAN 45 M SEC"* 
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Figure 4-5. PROBABILITY DISTRIBUTION CURVE OF STANDARDIZED WIND SPEC 

DIFFERENCE AVERAGED OVER ELAPSED TINES AT WAVE NUMBER 0.002 
CYCLES M-' FOR MIND SPEEDS LESS THAN 45 M SEC-1 
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PROBABILITY DISTRIBUTION CURVE OF STANDARDIZED WIND SPECTRA 
DIFFERENCE AVERAGED OVER ELAPSED TIMES AT WAVE NUMBER 0.0025 
CYCLES M-1 FOR WIND SPEEDS LESS THAN 45 M SEC-1 
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ngura 4-7. PROBABILITY DISTRIBUTION CURVE OF STANDARDIZED WIND SPECTRA 
DIFFERENCE AVERAGED OVER ELAPSED TIMES AT WAVE NUMBER 0.00375 
CYCLES M-l FOR WIND SPEEDS LESS THAN 45 M SEC* 1 
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Figure 4-9. PROBABILITY DISTRIBUTION CURVE OF STANDARDIZED WIND SPECTRA 
DIFFERENCE AVERAGED OVER ELAPSED TIMES AT WAVE, NUMBER 0.00625 
CYCLES M-1 FOR WIND SPEEDS LESS THAN 45 M SEC' 1 
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Figure 4-12. PROBABILITY DISTRIBUTION CURVE OF STANDARDIZED WIND SPECTRA 
DIFFERENCE AVERAGED OVER ELAPSED TIMES AT WAVE NUMBER 0.00025 
CYCLES M-l FOR WIND SPEEDS GREATER THAN 45 M SEC'l 


4-17 



I 



















Figure 4-15. PROBABILITY DISTRIBUTION CURVE OF STANDARDIZED WIND SPECTRA 
DIFFERENCE AVERAGED OVER ELAPSED TIMES AT WAVE NUMBER 0.0015 
CYCLES M-l FOR WIND SPEEDS GREATER THAN 45 M SEC-1 





NORTHROP SERVICES, INC. 


M-1259 



4-21 

















Figure 4-18. PROBABILITY DISTRIBUTION CURVE OF STANDARDIZED MIND SPECTRA 
DIFFERENCE AVERAGED OVER ELAPSED TIMES AT WAVE NUMBER 0.00375 
CYCLES M-l FOR WIND SPEEDS GREATER THAN 45 M SEC-1 





PROBABILITY (%) 

Figure 4-19. PROBABILITY DISTRIBUTION CURVE OF STANDARDIZED WIND SPECTRA 
DIFFERENCE AVERAGED OVER ELAPSED TIMES AT WAVE NUMBER 0.005 
CYCLES M-l FOR WIND SPEEDS GREATER THAN 45 M SEC"' 
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Figure 4-23. STANDARD DEVIATION OF WIND SPECTRA DIFFERENCE VERSUS ELAPSED 
TIMES AT VARIOUS WAVE NUMBERS FOR WIND SPEEDS LESS THAN 
45 M SEC-1 


4-28 
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Figure 4-28. STANDARD DEVIATION OF WIND SPECTRA DIFFERENCE AVERAGED OVER 
ELAPSED TIMES AND WIND SPEEDS VERSUS WAVE NUMBERS AND THE 

POWER LAW o (K) = o(K 0 )(£-) q WITH Kq = 0.001 CYCLES M' 1 AND 

q = -2.567 0 
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CONCLUSIONS 
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Figures 4-11 and 4-22 are the "universal" curves for the probability of the 
wind spectra difference. These two curves agree very well between one and 99 
percent level. The mean of the wind spectra difference 1 8 equal to zero. 

The standard deviation of the wind spectra difference can be found by using 
equation (3-1). By employing the universal curve, the wind spectra at later 
time can be estimated by giving present wind spectra. 

The universal curves show 57 percent of the time the spectral changes are 
negative. This indicates most of the time the wind spectra decrease. So the 
detailed wind spectrum on the average build-up faster than it decays in an 
Eulerian context. 



at wave number 0.001 cycles m ^ for wind speeds less than 45 msec From 

Figure 4-11, the 95 percent value of standardized wind spectra difference is 1.9. 

From Figure 4-28, the standard deviation of wind spectra difference for wave 

number 0.001 cycles m is 0,175. So the wind spectra difference at wave 

number 0.001 cycles m ^ is 0.3325, and the wind spectra at a later time will 

2 -1 -1 -1 

have 95 percent probability of not exceeding 10.3325 ra sec (cycles m ) at 
wave number 0.001 cycles m 
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